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In this work we study the consistency of a possible non-vanishing coupling HpLT of the order 
of 3.6 X 10~® as pointed recently by the CMS and ATLAS collaborations i. i , with measured 
lepton flavor violation processes involving quarkonium. We show that the most promising channel to 
confirm this excess is to look for the lepton flavor tan decay into a /o and p where the experimental 
limit could strongly improved with the new B factories as Belle II. 


Recently CMS collaboration has observed a slight excess of signal events with a significance of 2.4 standard deviations 
which can be interpreted as a Higgs particles decaying into a muon and tau leptons: 


Br{H — fir) 


0.84+°f7% 

0.77 ±0.62 % [2] 


( 1 ) 


Using the reported value of Higgs mass to be 125 GeV this requires the coupling Hfir to be of the order of 
3.6 X 10“^. Even if this observation is very challenging to be explained in new physics models [^, this value of the 
lepton flavor violating coupling is not in contradiction with the experimental upper limits for r —>■ /ry and r —^ 3/r. 
Indeed, using the formalism in Refs. [^, @ we obtain the values listed in Table [J for the corresponding branching 
ratios. The ejected value from H ^ ptr is done assuming that the Higgs couplings to charged leptons are given by 
SM values (3,0. 

The computation of Higgs-induced lepton flavor violation (LEV) in channel decays involving only charged leptons 
are even more intricated as the smallness of the lepton Yukawa couplings imply that higher loop contributions can 
be bigger than tree-level one [ 1 , 0 . To avoid this problem, we shall study the effect of the LEV H —>■ r/r coupling in 
processes involving quarkonium. In Ref. the effects of heavy sterile Majorana neutrinos in LEV decays of vector 
quarkonia has been studied. Sterile Majorana neutrinos induces 'jhlj, Zljj LEV couplings at one loop and WWUlj 
LEV couplings at tree level. These couplings produce LEV effects in quarkonia decay studied in Ref. 0- Another 
effect not studied there, is to produce a non-vanishing Hlilj LEV coupling at one level, whose effects in quarkonium 
decay were no analyzed in Ref. 0. 

In this paper, we assume no specific models for new physics behind the Hr/i coupling and we systematically study 
the LEV decays of quarkonia involving the Hth coupling, considering the phenomenological value of 3.6 x 10“^ 
pointed by the CMS and ATLAS collaborations. We shall show that even if the expected branching ratio are still 
below the experimental limit, some of them could be accessible to next generation of B factories as Belle H. 

The description of processes involving the annihilation or creation of heavy quarkonia can be systematically done in 
the framework of non-relativistic quantum chromodynamics (NRQCD) 0. This is a systematic expansion in terms 
of Qfs and the quarks relative velocity v with a clear separation of the perturbative phenomena occurring at the scale 
mq and the non-perturbative ones occurring at the scale mqv. The non-perturbative effects are encoded in universal 
matrix elements with a well defined hierarchy in the the v expansion. The novelty of this systematic approach is that, 
for some processes, color-octet configurations of the created or annihilated quark-antiquark pair yield contributions 
of the same order as the old color-singlet contributions to a given order in the as and v expansion. 

In this work we are interested in the order of magnitude of the branching ratios of the considered processes and 
will focus on the color singlet contributions which can be calculated using the old quarkonium techniques described 
in 00. A more refined analysis can be done in the most promising channels but this is beyond the scope of the 



Experimental bound [12,1 

Expected from 7/ —>■ pr 

r — >■ p7 

r —>■ 3p 

4.4 X 10“® 

2.1 X 10"® 

1.3 X 10“® 

1 X 10”“ 


TABLE I: LEV in r LEV decays involving charged leptons 
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present work. 

The invariant amplitude for the annihilation of color-singlet quarkonium in a Lj angular momentum configu¬ 
ration —>• X is given by MM 


M[QQ[^^+^Lj]^ X] = 


d'^q 

{2^Y 


Tr[0{Q,q)x{Q,q)\, 


where 0{Q,q) is the operator entering amplitude for the corresponding free quarks transition 
MlQi^ - q,S 2 ),Q{^ +q,si) -)• X] = v{^ - q,S2)0{Q,q)u{^ +g,si), 
and x{QjQ) denotes the wave function for the QQY^^^Lj] bound state 


( 2 ) 


(3) 


x(Q,g)= E - ^)YLM{ci)Ps,sAQ,q){LM■ssYJJz)■ 

M,s^ 


Here, Ps,s^ stands for the spin projectors 


Ps,sAQ,q) = \j^'^^ij + 9^^Av{j-q,s2){^si-^s2\ssA 


N, 


mq \ 2\/2mQ I 2 


i^ + ^ + mg) 


^2 


C^+i-mg) for 


5 ' = 0 
5=1 


(4) 


(5) 

( 6 ) 


where e(Q, SA) denotes the polarization vector of the spin one system. 

For s-wave quarkonium the wave function is rapidly damped in the relative momentum q and the leading terms are 
given by (Q, 0) and 0{Q, 0). In the zero-binding approximation the quarkonium mass M is given by M « 2mg 
and the amplitude reads 


R(0) 


vi[<3or«s.l^A'l = ^y—r. 

with M denoting the quarkonium physical mass and 

f (Aq 


{2^)- 


0{QA) 


■V'oo(q) = 


r 

m^sA 


m 

\fAA 


'-M) 


for 


5 = 0 
5 = 1 


(7) 


( 8 ) 


A similar calculation of the invariant amplitude for the production of color singlet quarkonium, X —>■ Sj] + Y 

yields 


i?(0) 


MiX^QQr-^SA+Y) = -^^-Tr 


OiQA) 


r 

AQ^sa 


-I- M) 


for 


5 = 0 
5 = 1 


(9) 


For p-wave quarkonium, the wave function at the origin vanishes and the leading term for the annihilation amplitude 
is given by the first term in the expansion in q of Eq. ©. A straightforward calculation yields 


M\QQf^+^Pj]^X] = -iJ2 AM-,SSYJJz)eAM)\l —R'iO)Tr [O^{Q,0)PsMQA) + OiQ,0)PAsAQA)] 

M,S^ 


where 


^Qct 


( 10 ) 

( 11 ) 


and in this case 


d^q 

(27r) 


3g“'*/'iM(q) = -i\l -^R'i0)eaiM). 


( 12 ) 
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The polarization vector £a{M) satisfies the following relations 


Y,{lM-lS,\m)eo.{M)ep{S,) = 

1 QaQ 13 

(13) 

M,S^ 


M,S^ 

{Jz)i 

(14) 

J2{^M;lS,\2J,)s^iM)£piS,) = 


(15) 


M,S^ 


The Higgs to Lj\ quarkonium coupling is obtained from the diagram in Fig. ([T]) which yields the following 

operator 


0{Q,q)=^^, 

V 


(16) 


where v stands for the Higgs vacuum expectation value. Using this operator in the previous formulae it is easy to 


Q , 

^+q,si 



&S+1L 


Q 


H 


% -q,s2 


FIG. 1: Feynman diagrams for the Higgs-quarkonium coupling. 

show that the only non-vanishing coupling of the Higgs to quarkonium is to S' = 1, J = 0 p-wave quarkonium, in 
which case we obtain 


M [QQfPo] ^H] = (17) 

V V TT 

Notice that this coupling is proportional to the derivative of the wave function at the origin, which according to the 
NRQCD rules is suppressed by a v'^ factor with respect to the wave function at the origin. This makes the radiative 
transitions involving s-wave quarkonium configurations of the same order as the non-radiative ones involving p-wave 
quarkonium configurations. The radiation changes the quarkonium quantum numbers allowing the corresponding 
quarkonium to couple to the Higgs. The calculation of Higgs-mediated lepton flavor violating radiative transitions 
involving s-wave quarkonium requires to work out Higgs-quarkonium-photon coupling. This transition is induced by 
the diagrams in Fig. ©■ 



FIG. 2: Feynman diagrams for the Higgs-quarkonium-photon coupling. 
From these diagrams, we identify the transition operator as 

f - ^ -I- ^ - mg 


0{Q,q) = ieeq^ 

V 


l + ^ + ji + ’ng 


{^ + q + k)^-ml (f-g-hfc)2-m^ 




where eg stands for the heavy quark charge in units of e. For s-wave J = 0 from Eq. ([7]) we obtain 


M[H ^QQ[^Soh] = -i 


.eeQMR{0) 


Av 


AttM 


Tr 


V(fc)(f+ i^)-(| + i^)^(fc)\ 


{% + kf-W?n 




-M) 


= 0 . 


(18) 


(19) 
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Similarly for s-wave J = 1 we get 


M[H ^QQ[-^Si]j]=i 


.eeQmQR{0) 


2v 


4:TTM 


Tr 


+ ^) — (^ + 

(f + fc)2-TO^ 




( 20 ) 


where r](Q) stands for the polarization vector of the quarkonium. A straightforward calculation yields 


M[H^ ggp5i]7] = J— 


( 21 ) 


with 


Tfiu — 9fiU 


Q^^k'' 
Q ■ k 


( 22 ) 


Now we focus on the Higgs mediated LFV processes. We start with the gQ[^Po] —>■ decay through the diagram 

in Fig ([3]). A direct calculation yields the following decay width 



FIG. 3: Diagram for the QQ^Po] —>■ decay. 


nQQfPo] ^ 9~T+] = 


27y^ \R'{0)\^M^ 
Stt^ v'^rnjj 



(23) 


where we neglected the muon mass and y stands for the Hyr coupling. 

Next we go through the corresponding radiative process QQpS'i] —)> /try. This decay proceeds through the diagrams 
in Fig. (gl). 



FIG. 4: Feynman diagrams for the —>■ yr decay. 


Neglecting the muon mass we obtain the following decay width 


F[ggp5i] ^ /r-T+y] 


ae^y^|i?(0)P \ 

127r3ti2 rnjj \M‘^) ’ 


(24) 


where 


f{x) = I — 6a: + 3a;^ + 2a:^ — 6 ln(a:). (25) 

The HyiT coupling can also mediate LFV decays of the tau meson involving light quarkonium. Although this is 
beyond the scope of the systematic NRQCD expansion due to the light quark mass, we still can use the quarkonium 
techniques taking care of extracting the corresponding non-perturbative pieces from the appropriate experimental 
data. The first possible decay is r —>■ /r gg[^Po] which goes through the diagram shown in Fig. ([5]). The decay width 
is given by 


272/21i?'(0)12 rrirM 



r(r- ^ y- ggpPo]) 


167r2'(;2 


(26) 
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FIG. 6: Feynman diagram for the r —>■ SiYl decay. 


where we neglected the muon mass. 

The corresponding radiative decay is r —?> p. QQ^Si\'y. The Feynman diagrams for this process are given in Fig. 
©• Neglecting the muon mass we obtain the following decay width 


where 


r(r- ^ QQfSi]-f) 


1-^(0)P mlM 

327r3u^ 



h{x) 


(1 + ^a;[(l 


a;)(3 — x) + 2 ln(a;)]. 


(27) 


(28) 


Finally, although the calculation does not require to use the quarkonium techniques, it is interesting to estimate the 
effects of the H^t coupling in LFV decay of gauge bosons. As a sample we calculate W —>■ The non-perturbative 
piece of this decay is related to the pion decay constant. This decay is induced by the diagram in Fig. ([7]). 



FIG. 7: Feynman diagram for the W —>■ r/in decay. 
The amplitude for W{Q,e) —>■ ^J.{pl)T{p 2 )^T{p 3 ) is 

.. yg'^VudU P3 -£{Q) ^ 

M = —y=- - ---r2- ^U{p2)v{pi), 

2y/2mw (Pi + P 2 ) - my 
where g is the weak coupling constant. The resulting decay width is 

2„2 \2t/2 


f2 2 V2y2 


(29) 


where 


/(a, b) = -^ [a^ {b^ - l) (24a® - 6a^ {4b^ + 7) + {b^ + 17) (26^ + l) - 66^) 

— 6 (a^ — l)^ (4a® — a"^ (66^ + l) + 2a^6^ + b'^) In — 126"^ln(5) 


( 30 ) 
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Process 

re:.p(GeF) 

|R(0)|^(GeF®) 

|R'(0)|^(GeF®) 

X ^ e+e“ 

1.28 X 10"® 

4.856 

- 

J/tp e'^e~ 

5.54 X 10"® 

0.560 

- 

4> —>■ e''"e“ 

1.26 X 10"® 

5.53 X 10"^ 

- 

xl 77 

2.34 X 10"® 

- 

3.10 X 10"^ 

/o ->■ 77 

0.29 X 10"® 

- 

1.08 X 10"“* 


TABLE II: Numerical values of the non-perturbative matrix elements extracted from the leptonic and two photon decays of 
quarkonia. 


with a = b = and we neglected the pion and muon masses. 

The results for the studied decays depend on the color-singlet matrix elements i?(0) for quarkonium and R'{0) 
for quarkonium. The same matrix elements appear in the leptonic decay of the first and two photon decays of 
the latter. We use the available experimental results on these decays to extract the phenomenological value of the 
matrix elements. The only matrix element that cannot be calculated this way is R'{0) for the ^Pq bb quarkonium Xbo- 
There is no available experimental data on the Xbo 77 transition, but R'{0) has been calculated in several potential 
models summarized in Ref. [l^ yielding R'{0) ~ 1 GeV^ and we will use this value in our calculations. 

The leptonic decay of vector quarkonia is induced by the diagram in Fig. ([5]). The corresponding decay width is 




% + q,si 

7 

1+ 

% -q,S 2 


r 


FIG. 8: Diagram for the QQ^Si] —>■ decay. 


r(QQp5i] ^ i+r) 


4a^el\Rm^ 

M2 


(31) 


where we neglected the lepton mass. The two photon decay of ^Pq quarkonium proceeds through the diagrams in 
Fig. ([S]). The decay width is 



FIG. 9: Feynman diagrams for the two photon decay of ^Po quarkonium. 


r(^Po 77) 


432a2e||i?'(0)|2 

M4 


(32) 


In Table HIl we give the results for the matrix elements of the different quarkonia. As a first approximation we use an 
ss configuration for the /o(980). There is no available information on the total width of the XbO: thus we report the 
decay width when we use ^^^,,(0) = 1 GeV^ from quark model calculations [l5l| : 


F(xm (J-t) = 5.5 X GeV. 


(33) 


The branching ratios of the remaining decays are calculated using the estimates for the non-perturbative matrix 
elements in Table HIl We list in Table HID the so obtained branching ratios. Here, the branching ratios include a factor 2 
to account for the two charge states where appropriate, e.g. BR{xcO M^) = BR{xcO —t +BR{XcO —t 

In general these branching ratios are small. The most promising decay is the r —>■ /i/p. We recall that we assumed 
an ss conhguration for this meson. The nature of the low lying scalar mesons is an old problem (see [l^ and references 
therein) and it would be desirable to have a closer approximation to the non-perturbative effects in this decay. 
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Process 

Branching Ratio 

Exp. bound 

Xco -5- /rr 

1.5 X 10"^'^ 


T —>■ firy 

5.7 X 10"^'* 


J/'4> — >■ /iry 

5.1 X 10"^'^ 


r ^/i/o(980) 

8.4 X 10"^^ 

< 3.4 X 10"® 

r —>■ /i()>7 

1.7 X lO"^"* 


W —>■ /iTTT 

3.2 X 10"^'^ 



TABLE III: Branching ratios for lepton flavor violation decays involving the Hfir coupling. 
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